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ABSTRACT 

We calculate Yukawa corrections of order aewMf /My^ to single top quark produc- 
tion via qq' — s- tb at the Fermilab Tevatron in the two-Higgs-doublet models. In 
our calculation we also keep the terms proportional to A/^ tan^ /3 since their effects 
may become rather important for large tan/3. The corrections can amount to more 
than a 15% reduction in the production cross section relative to the tree level result 
in the general two-Higgs-doublet model, and a 10% enhancement in the minimal 
supersymmetric model, which might be observable at a high-luminosity Tevatron. 
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1. Introduction 



The top-quark physics has become a very active research area since the top quark was 
discovered by the CDF and DO Collaborations at the Fermilab Tevatron[l]. It is timely to 
focus attention on directly investigating the properties of the top quark, especially its pro- 
duction mechanisms. With the increase in the number of top quark events at the Tevatron, 
the experimental errors are expected to be further reduced. With the next Tevatron run at 
•y/s = 2.0 TeV, one can expect about twenty times as much data as exist now. Thus, the 
comparison between the observed top quark production properties and more precise theoret- 
ical calculations will be an inportant probe for the possible existence of new physics. At the 
Tevatron top quarks are produced primarily via two independet mechanisms: The dominant 
production mechanism is the QCD pair production process qq — > tt[2\ . Single top production 
via 1^-gluon fusion subprocess g + W ^ tb [3] and the subprocess qq' — > tb [4] are also im- 
portant. These latter processes involve the electroweak interaction and, therefore, can probe 
the electroweak sector of the theory, in contrast to the QCD pair production mechanism. 
A recent analysis [5] of the process qq' — > tb showed that it is potentially observable at the 
Tevatron with 2-3 fb"^ of integrated luminosity. This process probes the top quark with a 
timelike W boson, > {Mt + M^y, while the 1^-gluon fusion process involves a spacehke W 
boson, q"^ < 0, and these processes are therefore complementary. Moreover, in the Standard 
Model (SM), the process qq' — > tb can be rehably predicted and the theoretical uncertainty 
in the cross section is only about a few percent due to QCD corrections [6]. Although the 
statistical error in the measured cross section for this process at the Tevatron will be about 
±30% [5], a high- luminosity Tevatron would allow a measurement of the cross section with a 
statistical uncertainty of about 6% [6]. At this level of experimental accuracy a calculation of 
the radiative corrections is necessary. In Ref . [6] the QCD and Yukawa corrections to single 
top quark production qq' — > tb have been calculated in the SM. While the QCD corrections 
were found to be quite large, the Yukawa corrections were found to be negligible. Since the 
SM weak corrections are expected to be comparable to the Yukawa corrections, they too 
should be negligible. Beyond the SM, the Yukawa corrections might be greatly enhanced, 
since more Higgs bosons with stronger couplings to top or bottom quarks are involed in 
some new physics models. Once the top quark mass is known precisely, these effects could 
be used as an indirect test for new physics beyond the SM; for example, the two-Higgs- 
doublet model(2HDM) and the minimal supersymmetric model(MSSM) [7] . At least, the 
data could be used to place restrictions on these models. Therefore, it is worthwhile to in- 
vestigate single top quark production via qq' tb in these models. In this paper we present 



the calculation of the Yukawa corrections of order a^.^M^ jM"^ to single top production at 
the Fermilab Tevatron in both the 2HDM and the MSSM. These corrections arise from the 
virtual effects of the third family (top and bottom) of quarks, neutral and charged Higgs 
bosons, and neutral and charged Goldstone bosons. We note that our calculations can be 
easily extended to the pseudo-Goldstone boson (PGB) corrections in technicolor models [8] 
by substituting the virtual PGB's in the technicolor models for the virtual Higgs bosons in 



The tree-level Feynman diagram for single top quark production via gg' — * ih is shown in 
Fig. 1(a). The Yukawa corrections of order a^.^M'^ / M'^r to the process qq' — » tb arise from the 
Feynman diagrams shown in Figs.lb-lh. In our calculations, we used dimensional regalariza- 
tion to control all the ultraviolet divergences in the virtual loop corrections and we adopted 
the on- mass -shell renormalization scheme [9]. We also kept all the terms proportional to the 
product Mfotan/? in the charged Higgs couplings to the third family of quarks since these 
effects may become rather important for large tan (3. We used the 't Hooft-Feynman gauge 
for the propagators of virtual W boson and Goldstone bosons. Including the 0{aewMf / M'^^) 
Yukawa corrections, the renormalized amplitude for qq' —>■ tb can be written as 



where Mq is the tree-level matrix element, and ^M^^*^'*^^ represents the 0{aeyjM^ /My^r) 
Yukawa corrections arising from the self-energy diagrams Figs, lb- le and vertex diagrams 
Figs. If- Ih. These are given by 



the 2HDM and MSSM. 



2. Calculations 



Mren = Mq + 5M' 



■vertex 



(1) 



Mo = i 




v{P2)'ynPLu{pi)u{p3)YPLv{p4) 




and 



5M' 




■v{p2h^.PLu{pi)u{p,)[rPLSi + rPR52 + AIPl + KPrHpa) (3) 



where 
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+ T.^'^l) finite + /l , 



(4) 
(5) 



P'3f2+P'ji, 
P'3f2^+Plf3^- 



(6) 
(7) 



Here pi and p2 denote the momentum of the incoming quarks q and q', while ^3 and ^4 are 
used for the outgoing t and b quarks, and s is the center-of-mass energy of the subprocess. 
SZj^ and SZ^ are the wave-function renormahzation constants, and //"'^ are form factors 
which are presented in Appendix A. 

The renormahzed differential cross section of the subprocess is 

da s — M?^. , , ,0 
dcosO 32ns^ ^ 

where Q is the angle between the top quark and incoming quark. Integrating this subprocess 
differential cross section over cos Q one finds 

a = ao + Aa (9) 

where 

= 4^ /{s - M^^Y ^^' - + - + ^^'^ + 

is the tree-level result and the correction is 

-[\{s - M^f + M^{s - M^) + Mt\M,f^ + \s{s - M^){M,fi - M,f^) 
+sMtMb{Mtf^ - Mbfi) + 2sMtMb62} . (11) 

The total hadronic cross section for the production of single-top-quark via qq' can be written 



in the form 



^(«) = E / dxidx2aij{xiX2S,M^,iJ,^)[ff{xi,ii)ff{x2,ii) + {A ^ B)], (12) 

where 

s = (P1 + P2)', (13) 

s = X1X2S, (14) 

Pi = XiPi, (15) 

and 

P2 = X2P2- (16) 



Here A and B denote the incident hadrons and Pi and P2 are their four-momenta, while 
i.j are the initial partons and Xi and X2 are their longitudinal momentum fractions. The 
functions ff^ and are the usual parton distributions. In our numerical calculations, we 
have used the CTEQ3L parton distribution functions for the tree level cross section, and 
CTEQ3M parton distribution functions[10] for the O^af^^M^ /My^,) Yukawa corrections, as 
in ref.[6], to facilitate comparison. There is no Yukawa correction to parton distribution 
functions as pointed out in Ref.[ll]. We will also compare our calculations with those 
calculated using the MRSG parton distribution functions[12] below. Finally, introducing the 
convenient variable r = X1X2, and changing independent variables, the total corss section 
becomes 



where tq = {Mt + Mb^/s. The quantity dLij/dr is the parton luminosity, which is defined 
to be 



In the following we present some numerical results for the Yukawa corrections to the total 
cross section for single top quark production via qq' tb at the Fermilab Tevatron with ^/s — 
2 TeV. In our numerical calculations we chose Mz — 91.188GeV, Mw = 80.33GeV, Mb — 
5GeV,aew — 1/128, and /i = VI. The Higgs masses Mh, Mh, Ma, Mh+ and parameters 
a, 13 are not constrained in the general Two-Higgs-Doublet Model, but in the Minimal SUSY 
Model, relations[13] among these parameters are required by supersymmetry, leaving only 
two parameters free; e.g., Mh+ and tan/?. Also, in the MSSM the charged Higgs mass is 
heavier than the W mass due to the relation = M,^ + Al\. The experimental lower 

limit for the charged Higgs mass is 44.1GeV[14], independent of the additional parameters 
a and j3. In our calculations we will use Mh+ = 600GeV for 2HDM and Mh+ > lOOGeV 
for MSSM. The upper bound of tan/?; viz, tan/i < 0.52GeV~^MH+, has been determined 
from data on B ^ rz/X[15]. The lower limits on tan/? are tan/? > 0.6 from perturbative 
bounds [16] and tan/3 > 0.25 (for Mt = 175GeV) from perturbative unitarity[16]. We will, 
therefore, limit the value of tan /? to be in the range of 0.25 to 30. 




(17) 




(18) 
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Numerical results and discussions 



Figure 2 shows the relative correction Act/cto as a function of M^j using the CTEQ3L 
parton distributions for the tree-level cross section (Tq and the CTEQ3M parton distributions 
for the correction Ac, as in ref. [6]. The solid curve corresponds to the 2HDM assuming 
Mh = Mh and Mt = 175GeV, Mh+ = Ma = GOOGeV and tan/5 = 0.25. The dotted curve 
corresponds to the SM for Mt = 175 GeV. For M^ = Mh the corrections are independent of 
a in the 2HDM. The corrections in the SM are negeligbly small, in agreement with Ref. [6]. 
However, in the 2HDM, the corrections can reduce the cross section by more than —10% for 
Mh = Mh < lOOGeV, and for Mh = Mh = 50GeV they can be as large as -20%. 

Figure 3 shows the relative correction Acx/ctq in the MSSM, assuming tan/3 = 0.25. Since 
the corrections are not sensitive to Mh+, for Mh+ > 400GeV, we only present the results 
for Mh+ in the range lOOGeV to 400 GeV. In Fig.3 the solid curve corresponds to Mt = 175 
GeV, again using CTEQ3L distributions for (Tq and CTEQ3M distributions for Acr. The 
dotted curve corresponds to Mf = 175 GeV but using the MRSG parton distributions. The 
dashed curve corresponds to Mt = 200 GeV and using CTEQ3L for ctq and CTEQ3M for 
Act. For a light charged Higgs, the corrections can be quite significant. For Mh+ = 100 
GeV the corrections reach 9% for Mt = 175 GeV and 13% for Mt = 200 GeV. The curve for 
Mt = 175GeV has a peak at Mh+ = 170 GeV and the curve for Mf = 200GeV has a peak 
at Mh+ = 195 GeV, which is due to the fact that M5 = 5 GeV and the threshold for open 
top decay into a charged Higgs plus a bottom is crossed in these regions. If we change the 
top quark mass, we found that this region is also shifted correspondingly, which provides 
a check on our calculations, especially of the treatment of the threshold. From Fig.3 we 
see that the difference between the results obtained using CTEQ3 distributions and using 
MRSG distributions is negeligibly small. We also found that the results using MRS (A') 
distributions [12] are almost the same as the MRSG results, and thus we did not present 
them. 

In both Fig. 2 and Fig.3, we used the minimal value (0.25) for tan/?. When tan/5 becomes 
larger, the corrections may drop rapidly since the dominant effects arise from the terms 
~ a^w ^*^2 p ■ In Fig. 4 we present the dependence of the relative correction, Ao"/o"o, on the 
value of tan/3 using CTEQ3L for (Tq and CTEQ3M for Aa. The solid curve corresponds to 
the 2HDM assuming Mt = 175GeV, Mh+ = Ma = 600GeV and Mh = Mh = 65GeV. The 
dotted curve corresponds to the MSSM assuming Mt = 175GeV and Mh+ = lOOGeV. The 
corrections are only significant for small tan f3 and are quite sensitive to tan f3 for tan f3 < 5. 

Since the cross section for single top production can be reliably predicted in the SM 
[6] and the statistical error in the measurment of the cross section will be about 6% at a 



high-luminosity Tevatron[6], these corrections may be observable; at least, interesting new 
constraints on these models can be established. 

Note that in the MSSM, besides these Yukawa corrections arising from the Higgs sector, 
the supersymmetric (SUSY) corrections due to super particles (sparticles) should also be 
taken into account [17]. The dominant virtual effects of sparticles arise from supersymmet- 
ric QCD corrections of order and the supersymmetric electroweak correction of order 
UewM^ / M^r which arise from loops of charginos and neutralinos, the supersymmetry part- 
ners of Higgs and vector bosons. However, the anomalous magnetic moment for a spin 1/2 
fermion vanishes in the SUSY limit [18] and away from the SUSY limit the cancellations 
have somewhat less effect. Therefore, in general one can expect the Yukawa corrections from 
the Higgs sector and the supersymmetric electroweak corrections from virtual charginos and 
neutralinos to cancel to some extent. 

In conclusion, we calculated the Yukawa corrections of order a^wM^ / M^^ to single top 
quark production via qq' tb at the Fermilab Tevatron in the general two-Higgs-doublet 
model and the minimal supersymmetric model. For favorable parameter values, the correc- 
tions can reduce the cross section by more than 15% in the general two-Higgs-doublet model 
and enhance the cross section by up to 10% in the minimal supersymmetric model. These 
effects may be observable at a high-luminosity Tevatron. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant No. DE-FG02-91-ER4086. 

Appendix A 



5Z^ 



i=A,G'> 



5Z^ 




A 



(bbi) 



i=A,GO ^ 
i=H+,G+ 



where A = ^ — 7£;-|-log47r with 7^ being the Euler constant and D = 4 — 2e is 
dimension. The functions F^^^\ G^^^^^ are defined by 



Fi^^'^) = Cdyy^ log 
Jo 



Mfy{y-l) + M]{l-y) + Mly 



y 



"+'(1-1/) 



Mfy{y-l) + Mm-y) + Mly 



(A.2) 
the space-time 

(A.3) 
(A.4) 



The form factors /^^ are obtained by 



L,R 



f,^j'^ are given by 



/21 

J 31 

hi 

fR 
J21 

fR 
J 31 



fl2 



f'. 
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/32 

fR 

J 12 

fR 
J 22 



MtM^^ijXj[rii{2c23 + C12 - 2C22) + Ai(-ci2 + 2c22 + 2c2i + cn 
MtM^^ijXj[r]i{-Ci2 - 2C22) + \{-Cn + 2c22 - 2c23 + C12)] 

2MtMb^ijXjr]iC24 
MSjXj[M^\{2c23 + C12 - 2C22) 

+M^7]i{-Ci2 + 2C22 + 2C21 + Cii - 4C23)] 

M6ei,-A,[M2Ai(-Ci2 - 2C22) + M^rii{-cu + 2c22 - 2c23 + C12)] 
2M^^ijXjXiC24: 

-MtM^^ijXj[r]i{2c23 + C12 - 2C22) 

-Ai(-Ci2 + 2C22 + 2C2l + Cii - 4C23)] 

-MtM^^,jXj[r]i{~ci2 - 2C22) - Xi{-cn + 2c22 - 2c23 + C12)] 
-2MtMS,\jV^C24 
Mb^ijXj[M^Xi{2c23 + C12 - 2C22) 
-M^r]i{-ci2 + 2c22 + 2c2i + Cn - 4C23)] 



(A.5) 



(A.6) 
4c23)] (A.7) 
(A.8) 
(A.9) 

(A.IO) 
(A.ll) 
(A.12) 

(A.13) 
(A. 14) 
(A.15) 

(A.16) 



fs2 = M,e.,A,- [M^\{-Ci2 - 2c22) - M^r]i{-cu + 2c22 - 2c23 + C12)] (A.17) 

/i^3 = -2M2ei,-7?,r?,C24 (A.18) 

/2^3 = -MtCijr]j[M^r]i{2C21 + 2c22-4:C23 + Ci2-Cu-Co) 

+M^Xi{2c23 - 2c22 + C12)] (A.19) 

/3^3 = -MSjr]j[M^r],{2c22-2c23 + 3ci2-Cn + Co)-M^X,{2c22 + Ci2)] (A.20) 

/i«3 = -2M^M^,e^,r/,AiC24 (A.21) 

= -M2Mb6,r/,[r/i(2c23-2c22 + Ci2) 

+Ai(2c2i + 2c22 - 4c23 + C12 - cii - Co)] (A.22) 

/3^ = -M^MhCijr]j[-r]i{2c22 + Ci2) + Xi{2c22-2c23 + 3ci2-Cu + Co)] (A.23) 

/i^4 = 2M^^ijr]jr]iC2A (A.24) 

/2^4 = -MSjrij{M^ri,{Ac23 - 2c2i - 2c22 + 3ci2 - 3cn - Cq) 

+M2A,(2c23 - 2c22 + C12)] (A.25) 

/si = -Mt^^jVj [Mhi{2c23 - 2c22 + C12 + Cn + Co) - M2Ai(2c22 + C12)] (A.26) 

/i^ = -2MtM6eyr?jAiC24 (A.27) 

/2^ = -M,2M6ei,-77,[?7i(2c22-2c23-Ci2) 

+Ai(2c2i + 2c22 - 4c23 + 3cii - 3ci2 + Co)] (A.28) 

/3^ = -M^Mb^.,jr]j \ni{2c22 + C12) + Ai(2c22 - 2c23 - C12 - Cn - co)] (A.29) 

= -M,2M,\A,(co + 2ci2-Cn) (A.30) 

/2^5 = 2MtM,\A,(ci2 + C23) (A.31) 

/3^5 = 2MtM,\A.(c,2 + c22) (A.32) 

/l^ = MtMb?7jAj[2c24 + s(ci2 + C23) -Mj2(co-C2l + Ci2 + C23)] (A.33) 

/2^5 = 2M^Mt,r)jXjico-C2i) (A.34) 

/3^5 = 2M2M,r7,A,(ci2 + C23) (A.35) 

ft, = M^M^{co + Cu) (A.36) 

/2^6 = 2M4M2(ci2 + C23) (A.37) 

/3^6 = 2MX(C12 + C22) (A.38) 

/l^ = MtM6[2c24 + s(ci2 + C23)-Mt2(ci2 + C23-C2l-2cii-Co)] (A.39) 

= -2M,2M6(co + 2cn + C2i) (A.40) 

ff, = 2M^M,icu + C23) (A.41) 



(A.42) 

where the sums over i = H+, G+, j = H,h for /^'^ and /^'^, i = H+, G+, j = A, for /^'^ 
and /i4'^, and j = H,h for /j^'^ are imphed. The functions defined as 

Cij = Cij{-ps,k,M,,M,,Mj) for /^'^ /^'^ 

% = Q,(-p3,A;,M,,M,-,M,) for 

= Cij{-ps,-p4,Mt,Mj,Mh) for /^'^^ 

Cij = Cij{-p3,-p4,Mt,MA,Mb) for /^^'-^ 

are the three-point Feynman integrals [19]. The constants rji, Aj and are defined as 



r]H+ ^ Va^ cot/3, 77G+ = rjGo = 1, (A.43) 

TjH — sin a/ sin /3, 77/1 = cos a/ sin /3, (A. 44) 

Aif+ = Aa = tan/3, Ag+ = Ago = 0, (A.45) 

\h = cos a/ cos /3, A/j = — since/cos/?, (A. 46) 

^H+H = -^G+h = sm{p-a), (A.47) 

^H+h = ^G+H = -cos{P - a), (A.48) 

Ch+a = Cg+go = 1, C//+GO = Cg+a = (A.49) 
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Figure Captions 



Fig.l Feynman diagrams: (a) tree-level, (b)-(e) self-energies, (f)-(li) vertex corrections. 

Fig. 2 The relative correction Aa/o-Q as a function of Mh using the CTEQ3L parton dis- 
tributions for o"o and the CTEQ3M parton distributions for Ao". The solid curve corresponds 
to the 2HDM assuming My, = Mh, Mt = 175 GeV, Mh+ = Ma = 600GeV and tan/? = 0.25. 
The dotted curve corresponds to the SM for Mj = 175 GeV. 

Fig. 3 The relative correction Aa/ao as a function of Mh+ in the MSSM, assuming tan (3 = 
0.25. The solid curve corresponds to Mt = 175 GeV using the CTEQ3L parton distributions 
for (To and the CTEQ3M parton distributions for Ao". The dotted curve corresponds to 
Mt = 175 GeV using the MRSG parton distributions. The dashed curve corresponds to 
Mt = 200 GeV using the CTEQ3L distributions for ctq and the CTEQ3M distributions for 
Act. 

Fig. 4 The relative correction Act/cto as a function of tan/3 using the CTEQ3L distribu- 
tions for (To and the CTEQ3M distributions for Aa. The solid curve corresponds to the 
2HDM assuming Mt = 175GeV, Mh+ = Ma = 600GeV and Mh = Mh = 65GeV. The 
dotted curve corresponds to the MSSM assuming Mt = 175 GeV and Mh+ = 100 GeV. 



